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Abstract 

We study a generalization of the "shear-free part" of the Goldberg-Sachs theorem for Einstein spacetimcs 
admitting a non-twisting multiple Weyl Aligned Null Direction (WAND) £ in n > 6 spacetime dimensions. 
The form of the corresponding optical matrix p is restricted by the algebraically special property in terms 
of the degeneracy of its eigenvalues. In particular, there necessarily exists at least one multiple eigenvalue 
and further constraints arise in various special cases. For example, when p is non-degenerate and the Weyl 
$H . components <&ij are non-zero, all eigenvalues of p coincide and such spacetimes thus correspond to the 



Robinson- Trautman (RT) class. On the other hand, in certain degenerate cases all non-zero eigenvalues can 
be distinct. We also present explicit examples of Einstein spacetimes admitting some of the permitted forms 
of p, including examples violating the "optical constraint" . The obtained restrictions on p are, however, 
in general not sufficient for £ to be a multiple WAND, as demonstrated by a few "counterexamples". We 
also discuss the geometrical meaning of these restrictions in terms of integrability properties of certain null 
distributions. Finally, we specialize our analysis to the six-dimensional case, where all the permitted forms of 
p are given in terms of just two parameters. In the appendices some examples are given and certain results 
pertaining to (possibly) twisting mWANDs of Einstein spacetimes are presented. 
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1 Introduction 

In four-dimensional General Relativity, a fundamental connection between geometric optics and the algebraic 
structure of the Weyl tensor is provided by the Goldberg-Sachs (GS) theorem [1,2], which states that an Einstein 
spacetime is algebraically special if, and only if, it contains a shearfree geodetic null congruence (cf. [3, 4] for 
related results and generalizations). This theorem plays an important role in the construction of algebraically 
special exact solutions of the Einstein equations [3], as the remarkable discovery of the Kerr metric shows [5]. 

In recent years, the growing interest in higher dimensional gravity has motivated the study of extensions 
of the above concepts to n > 4 dimensions. An algebraic classification of the Weyl tensor based on the notion 
of Weyl Aligned Null Directions (WANDs) has been put forward in [6] (see also [7] for a recent review). Fur- 
thermore, a higher dimensional version of Newman-Penrose (NP) and Geroch-Held-Penrose (GHP) formalisms 
have been presented in [8-10] and [11], respectively. However, simple examples reveal that neither the "geodesic 
part" nor the "shearfree part" of the GS theorem extend in an obvious way to higher dimensions (see [12] for 
a brief summary, along with the original references [8, 10, 13-17]). 

The proper formulation of the geodesic part of the higher dimensional GS theorem has been proven in [18] 
(see also [15, 17]): in particular, an Einstein spacetime admits a multiple WAND if, and only if, it admits 
a geodesic multiple WAND - hence one can restrict to geodesic multiple WANDs (mWANDs) with no loss 
of generality. Very recent work [12] has analyzed the shearfree part in five dimensions by proving necessary 
conditions on the form of the optical matrix p (defined below in (1.2)) following from the existence of an 
mWAND. 1 Contrary to the 4d case, the conditions obtained in [12] are, in general, not sufficient. In fact, it 
seems that in higher dimensions conditions that are both necessary and sufficient in general do not exist. 

In the present contribution we focus on n-dimensional Einstein (including Ricci-flat) spacetimes of type II 
(type D being understood as a special subcase thereof) in more than five dimensions and work out the corre- 
sponding necessary conditions on p under the assumption that the mWAND is twistfree (i.e., p is symmetric, 
equivalent to A; L j = 0, cf. (1.4) below). 2 Since the types N and III were already studied in [8], we will simply 
connect the results of [8] to our analysis when appropriate. As it turns out, in n > 5 dimensions the constraints 
on p are in general not as strong as those for n = 5 [12]. Nevertheless, combining our results for a type II 
spacetime with the results for type III and type N [8] gives 

1 See [19] for a different formulation of the "shearfree" condition and [12] for a comparison between the two approaches. 
2 This automatically guarantees that £ is geodesic. Except for a few remarks, we assume that p is non-zero, since we are interested 
in studying its possible non-trivial forms. 
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Theorem 1.1 (Eigenvalue structure of p for n > 6 and = 0). In an algebraically special Einstein spacetime 
of dimension n > 6 £/ia£ is not conformally flat, the (symmetric) optical matrix p of a non-twisting multiple 
WAND has at least one double eigenvalue. In the following special cases stronger conditions hold and the most 
general permitted eigenvalue structures are 

(i) /M-,] / 0: {o,o,0,...,0} [12], 

(ii) if det p ^ 0, &ij 7^ 0: {a, a, . . . ,a} (Robins on- Trautman, <&ij oc 5ij, type D(bd)), 
(Hi) if®ij = (type Il(abd)): {a, a, b, b, a, . . . , c n _ 6 }, 

(iv) if the type is N or III: {a, a, 0, . . . , 0} [8]. 

(The Weyl components <J>4 and &ij are defined below in (1.7).) The proof of (ii) and (iii) is a new result 
and will be given in sections 4 and 5, respectively. Note that in the above cases (i)-(iv) the matrix p possesses 
at least two double eigenvalues. In particular, in six dimensions all possible cases can be listed explicitly: 

Corollary 1.2 (Eigenvalue structure of p for n = 6 and Aij = 0). In a 6d algebraically special Einstein 
spacetime that is not confomally flat, the (symmetric) optical matrix p of a non-twisting multiple WAND can 
have only one of the following eigenvalue structures (where a,b might coincide, or vanish): (i) {a,a,b,b}; (ii) 
{a, a, b, 0}; (iii) {a, b, 0, 0}. If the spacetime is type III or N then the structure is {a, a, 0, 0}. 

Note that we do not claim that all classes compatible with theorem 1.1 or corollary 1.2 are non-empty. 

The structure of the paper is the following. In section 2 we briefly summarize results for general (possibly 
twisting) type II spacetimes in arbitrary dimension (already contained in [12]), in particular the so called 
"optical constraint". These will be useful in the rest of the paper, where, however, we limit ourselves to the 
non-twisting case. In section 3 we set up the study of general non-twisting type II spacetimes. It turns out 
that it is convenient to study the case of non-degenerate and degenerate p separately. This is done in sections 4 
and 5, where we derive certain necessary conditions for p being compatible with a non-twisting mWAND, in 
particular proving thus theorem 1.1. Using various explicit examples we show in section 6 that the necessary 
conditions obtained in sections 4 and 5 are in general not sufficient. In section 7 we elucidate the geometrical 
implications of theorem 1.1 in terms of existence of integrable null distributions with totally geodesic integral 
null two-surfaces. In section 8 we study the six-dimensional case (corollary 1.2) and provide various explicit 
examples of such metrics. 

Some related results are given in the appendices. In appendix A we present new algebraic constraints for 
general type II Einstein spacetimes. Appendix B contains a new proof of the canonical form of the optical 
matrix p for type III Einstein spacetimes in the case of a non-twisting mWAND £ (the original, longer proof 
was given in [8] ) . In appendix C we summarize various equivalent formulations of the "geodesic and shearfree" 
condition in four dimensions, and we discuss some extensions of those to higher dimensions (including optical 
structure, optical constraint, integrability of certain null distributions). Appendix D contains the proof (mostly 
relying on results of [12,20]) of the existence of an optical structure in all algebraically special spacetimes in 
five dimensions. In appendix E we present examples of twisting but shearfree mWANDs (in six dimensions). 
Finally, in appendix F examples of Einstein spacetimes of type D with non-degenerate mWANDs violating the 
optical constraint are given. 

Notation Throughout the paper we use higher dimensional GHP formalism developed in [11]. We employ a 
null frame 

{£ = e (0) = e« n = e (1) = e^,m {l) = e (i) = e«} (1.1) 

with indices i,j,k,... running from 2 to n — 1. The null vector fields £ and n and the orthonormal spacelike 
vector fields obey £ n = 1, • m^j) = and I • = = n • m^y The optical matrix p is defined 
by 

Pij = m a (i) m b {j) V b la, (1.2) 
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and its trace gives the expansion scalar (up to normalization) 

p = pu- (1.3) 

Its (anti-)symmetric parts are, respectively, 

Sij = P(ij), Aij = p^]. (1.4) 

We also define the rank of p as 

m = rank(p). (1-5) 
Other Ricci rotation coefficients used in the paper are [8, 10, 11] 

Lu = n a m b ^V b £ a , n = L a = mfyn b V b £ a , Mjo = m^ j} £ b V b m^ i)a , M jk = m^m^Vftm^. (1-6) 
Boost weight (b.w.) zero components of the Weyl tensor 

®ijki = Cijki, &ij = Conj, 2®fj = C iij, = Coioi (1-7) 

are subject to 

®ijkl = ®[ij][kl] = ®klij, ®i\jkl] = °' (I- 8 ) 

&ij = = ~\®ikjk, &tj = %]> * = ( L9 ) 

and thus they are determined by the ^(n — l)(n — 2) 2 (n — 3) components of &ij k i and the |(n — 2)(n — 3) 
components of <&fj. In a type II spacetime, at least one of the latter two quantities must be non- vanishing 
when £ is an mWAND (otherwise the spacetime would be type III or more special). Recall that the algebraic 
subtypes 11(a), 11(b), 11(c) and 11(d) are defined by the vanishing of <£, the symmetric traceless part of the 
"Weyl part" of &ijki, and <3?^, respectively (see [7]). 

There is no summation in case one or both repeated indices is/are in round brackets, unless it is said 
otherwise. E.g., there is no summation in &i(j)k(j) an d there is summation in <&ijkj and Ylj Q^^y 

In what follows we shall employ the Sachs equation [8, 10, 11] 

\>Pij = -PikPkj, (1-10) 
and the following Bianchi identities (eqs. (A10), (All) of [11]) 

p^- = -($ ifc + 2$f fc + $S ik )p kj , (1.11) 

-\><5> ijk l = 4®tjP[kl] ~ 2®[k\iPj\l] + 2$[k\jPi\l] + 2*ij[k\mPm\l]- ( L12 ) 

In a parallelly propagated frame, the GHP directional derivative b reduces to the NP derivative along the I 
direction, D = £ a V a . For an affinely parametrized I with an affine parameter r this is simply D = d r . 

There are two possible approaches to studying consequences of the Bianchi and Ricci equations. The first 
approach consists in applying the derivative operator p on certain algebraic equations and in deriving new 
algebraic constraints, e.g. eq. (2.3), using (1.10) (1.12) (cf. also appendix A and [12] (section 3.1 and appendix 
A therein)). In this paper we mainly use the second approach, which consists in solving the Sachs and Bianchi 
differential equations (1.10)-(1.12) and then in analyzing the compatibility of their solutions with the algebraic 
equations (2.1)-(2.3). 
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2 Results for general type II 



This section is devoted to summarizing some useful results that hold for all type II Einstein spacetimes, without 
assuming that p is non-twisting. In particular we present the algebraic restrictions on the Weyl components of 
b.w. zero that follow from the Bianchi identities, to be used in the next sections. More general new results are 
given in appendix A for future reference. For genuine type II spacetimes the (unique) mWAND is necessarily 
geodesic, while for type D spacetimes there always exists a geodesic mWAND [18]. In the frame used below 
we thus take I to be a geodesic mWAND, without loss of generality. 

2.1 Algebraic constraints 

For type II Einstein spacetimes, the Weyl tensor must obey the algebraic constraints (A. 12) and (B8) of [11] 
(already discussed in [15,17]) and b(B8) (derived in [12]), namely 



2 ®[jk\Pi\i] ~ 2 ®i{jPki] + &im\jk\Pm\q = 0, (2.1) 

$kjPij - <S> jk p i:j + QijPkj ~ <&jipjk + 2$ijpjk ~ <&ikp + <$>pik + QijklPjl = 0, (2.2) 

(2$fci - $jfc) Pupji + (2&ij - <5>ji) pjipkt - <5>%pjiPji + QpuPki + ®ijklPjsPis = 0. (2.3) 

Eq. (2.2) is traceless and its symmetric and antisymmetric parts read, respectively, 

(2$ fci - $ jfc ) Sij + (2$ - - S jk - $%p + <S>S ik + ^ijkiSji = 0, (2.4) 

QjkAji + QjiAkj + <5>ijp jk - $kjpji + <5>kiP + $A ik + ^ijuAji = 0. (2.5) 



The reason for investigating algebraic conditions such as (2.1)-(2.5) (or those of [8] in the case of type III/N) 
is that they will in general constrain the possible form of p. This way one obtains the standard Goldberg- Sachs 
theorem in four dimensions and one can arrive at similar conclusions also in higher dimensions, at least with 
additional assumptions (e.g. on the Weyl type [8], on the number of dimensions [12], on the form of the line- 
element [21,22], on the asymptotic behaviour [23], and/or on optical properties of £, as we shall discuss in the 
next sections). 

2.2 The optical constraint 

It was observed in [12] that the above conditions on the Weyl tensor appear to be less stringent when p satisfies 
the optical constraint [21] 

PikPjk^ P(ij)- (2-6) 

In particular, when this holds eq. (2.3) is not an extra restriction. This thus suggests that the branch of type 
II solutions whose mWAND obeys (2.6) corresponds to the case of a "generic" Weyl tensor [12]. It has been 
proven that (2.6) indeed holds for the Kerr-Schild vector of all (generalized) Kerr-Schild spacetimes [21,22] 
(including Myers-Perry black holes [13]), for non-degenerate geodesic double WANDs in asymptotically flat 
type II spacetimes [23] (see eq. (14) therein) and for the unique double WAND of all genuine type II Einstein 
spacetimes in five dimensions (see footnote 3 for the type D case) [12,20]. Nevertheless, double WANDs violating 
the optical constraint also exist, as shown by explicit examples constructed below in section 5.6.1 for n > 6, in 
appendix F for n > 7, and in section 6.3 of [12] for n = 5. 3 Although the above analysis applies only to type 
D and genuine type II spacetimes (eqs. (2.1)-(2.3) become trivial identities for more special types), it is worth 

3 We observe that all such examples are of type D. In fact, in five dimensions all type D Einstein spacetimes admitting a 
geodesic mWAND violating the optical constraint are known [12] (and coincide with the class of Einstein spacetimes admitting 
a non-geodesic mWAND [18]). Such "exceptional" null directions are necessarily twisting in 5d (but not in higher dimensions). 
However, in all such 5d type D spacetimes there always also exists a pair of non-aligned (non-twisting) mWANDs that do obey the 
optical constraint [12]. On the other hand, this is generically not true when n > 5 (see appendix F.5 for a ten-dimensional example 
admitting exactly two mWANDs, both violating the optical constraint). 
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remarking that the optical constraint is also obeyed by the mWAND of all type N Einstein spacetimes [8, 11]. 
The mWAND t of type III Einstein spacetimes also obeys the optical constraint provided either [8]: (i) the 
spacetime is five-dimensional; (ii) the Weyl tensor satisfies a genericity condition (see [8]); (iii) t is non-twisting 
(see also appendix B for a simpler proof in case (iii)). 

The optical constraint implies that (1 — ^p) is an orthogonal matrix [12] and that consequently [p, p T ] 
vanishes [12,21]. The optical matrix p is therefore a normal matrix and can thus be put, using spins, into a 
convenient block-diagonal form (see [21,24] for extended related discussions), i.e., 

,0,...,o). (2.7) 

The block-diagonal form (2.7) is useful for practical purposes because it allows for determining the r-dependence 
of p by integrating the Sachs equation (1.10) (in a parallelly transported frame) [24]. In the non-twisting case 
it reduces to a sequence of Is followed by a sequence of 0s (up to an overall factor). Note that the symmetric 
part of each two-block is proportional to a two-dimensional identity matrix, i.e. it is "shear-free". In four 
dimensions the optical constraint implies that either there is a single such block (in which case pij is shearfree) 
or that pij is symmetric with exactly one non-vanishing eigenvalue. However, the Goldberg-Sachs theorem 
shows that the latter case cannot occur. Therefore in 4d the optical constraint is a necessary condition for £ 
to be a repeated principal null direction but it is not sufficient [12]. 

On the other hand, in higher dimensions a vector field I obeying the optical constraint is in general shearing, 
as in the case of Myers- Perry black holes [13] (see also [14, 15] for a discussion of their optical properties). 
Together with other results [8, 10, 25], this has made clear that the shearfree condition is in general "too 
restrictive" in higher dimensions, as opposed to the four-dimensional case. In particular, it was observed in [10] 
that in odd dimensions a twisting geodesic mWAND is necessarily shearing. By contrast, twisting geodesic 
mWANDs with zero shear are permitted in even dimensions and they have necessarily det(p) / (as can be 
easily seen in a frame adapted to Aij, using the fact that Sij oc 5ij). An explicit example in six dimensions is 
discussed in appendix E. 

In the non-twisting case, the existence of shearfree spacetimes has been already known for some time in all 
dimensions - they are either RT [25] or Kundt [26] solutions, according to the presence/absence of expansion. 

2.3 Possible generalizations of the geodesic&shearfree property 

In arbitrary dimensions various geometric conditions can be considered which are different from the standard 
geodesic&shearfree condition (considered above) and which, however, become all equivalent (except for the 
optical constraint, cf. [12]) in the special case n = 4. Further such conditions are discussed in appendix C. One 
could thus conceive that various formulations of the Goldberg-Sachs theorem are in principle possible when 
n > 4. Some of these formulations have been studied in [12,19,27] (see also [21]) but none of these gives 
necessary and sufficient conditions. In the rest of this paper we will discuss necessary conditions determined 
by the presence of a non-twisting mWAND in an ra-dimensional Einstein spacetime and we will present several 
explicit examples. A possible interpretation of these results in terms of the geometric conditions of appendix C 
will be discussed in sections 7 (for n > 6) and 8.3 (for n = 6). 

3 Non-twisting £: general properties 

Here we study the optics of a hypersurface orthogonal mWAND I in a type II Einstein spacetime. This 
mWAND is automatically geodetic. Thus we have 

Ki = = A ij , (3.1) 

so that p^ = Sij. Since the algebraic equations (2.1)-(2.5) are trivial for Kundt spacetimes, in what follows 
we assume p 



p = adiag l,...l, 



1 + a 2 b 2 



1 — ab\ 1 1 —ab v 

ah 1 J ' ' ' ' ' 1 + a 2 b 2 [ab u 1 
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3.1 Case $A ^ 

This case has been already analyzed in [12], arriving at point (i) of theorem 1.1 of section 1. As already 
observed there, in this case p satisfies the optical constraint and it is necessarily degenerate (m = 2). Moreover 
it is shearing, except for p = (which is necessarily the case for n = 4). 
It remains to consider the case when <E>^ = 0. 

3.2 Case $A = 

This case defines the subtype 11(d) in the notation of [6]. For <E>^ = 0, eq. (2.2) reduces to (recall Ay = here) 

-p^>ik + <S>pik + ^ijpjk + ^ijkipji = 0, (3.2) 

and eq. (2.5) to 

[p,* 5 ]=0. (3.3) 

Note that taking into account (3.3), eq. (3.2) is symmetric and corresponds to eq. (2.4). 
Similarly, eq. (2.3) yields 

-QikPjlPjl + $(p 2 )ik + 2$ij(p 2 )jk + $ijki(p 2 )ji = 0. (3.4) 

Thanks to (3.3), one can choose a basis where both pij and §ij are diagonal, pij = diag(/92, P3, ■ ■ ■ ), 
$>ij = diag(3>2, $3, • • • )j therefore <&ijkj = for i ^ k. In this frame, the off-diagonal components of the 
algebraic constraint (3.2) are 

2 P(.3)®i(j)kU) = ( fOT k ^ *)■ ( 3 - 5 ) 

j 

The diagonal part of (3.2) can be expressed as Cijpj = 0, where pi is the vector pi = (p2,p3 • • • ) and the linear 
operator Cij is given by 

$2 $2 ••• $2 
$3 $3 ... $ 3 

Cij = Wij + + 2diag($ 2 , $3, • • • $n) 

$ n 

or, equivalently, £jj = Wjj + ($ + 2$(^)5jj — <3?(j), where we have defined the symmetric and traceless matrix 

= C(i)(j)(i)(j) ( no summation). (3-6) 

Note that W^) = 0. From the first of (1.9), it follows 

^Wi, = -2$i, (3.7) 
i 

which will be used in certain calculations throughout the paper. (From (3.7) it follows that Wjj = => Cij = 0.) 
Thus the sum of the rows of Cij vanishes and so they are linearly dependent and 

det£ = 0. (3.8) 

Therefore zero is an eigenvalue of C and (the diagonal part of) (3.2) admits non-trivial solutions. 

Note that (3.4) has the form Cijfj = 0, where components of the (n — 2) dimensional vector p 2 are squares 
of components of p. The characteristic polynomial of C is 

det(£ - XI) = \ n ~ 2 + £v t _ 3 A n - 3 + • • • + foA + k , (3.9) 
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with Uq = due to (3.8). Now zero is a multiple eigenvalue of C iff k\ = 0. Let us observe that for a generic 
form of a type II Weyl tensor k\ is non-vanishing and thus zero is a single eigenvalue of £. 4 If this is the case 
then p 2 is proportional to p and therefore 

pij = adiag(l, 1, ... 1, 0, ... 0). (3.10) 

Recalling also point (i) of theorem 1.1, we can conclude with 

Proposition 3.1. For non-twisting type II Einstein spacetimes with k\ / 0, all non-vanishing eigenvalues of 
the optical matrix p^ coincide. 

If <&(j / 0, the stronger result of point (i) of theorem 1.1 holds. 

As we will see in sections 5.6 and 8.1 (table 2), there do exist non-twisting type II Einstein spacetimes 
with distinct non-vanishing eigenvalues of p. According to Proposition 3.1, for these spacetimes k± vanishes. 
In particular C, and therefore also k\, vanishes for spacetimes with Wij = 0. We will study special cases with 
various parts of the Weyl tensor vanishing in the following sections. 



4 Non-twisting £: non-degenerate p (m = n — 2) 

First, let us consider the case of a non-degenerate p (i.e., detp 7= - this is relevant, e.g., for asymptotically 
flat algebraically special spacetimes [23]). When detp / one necessarily has (see point (i) of theorem 1.1) 

$tj = 0, (4.1) 

so that the type is 11(d). As noticed above we can use a frame in which both p^ and <E>jj are diagonal, which 
is moreover compatible with parallel transport [24,28]. Thus the eigenvalues of the optical matrix are [28] 

In order to discuss the possible structures of p it is convenient to discuss separately various cases in which 
different parts of the Weyl tensor are (non-)zero. 

4.1 Cases $^ 7^ and = 7^ 

From Bianchi equation (1.11) one has 

d(^-\ =-<$>. (4.3) 

V P(i) J 

The case when all pi coincide (and are non-zero) is the known RT case, for which also all <£j must coincide [25,28]. 
Thus let us consider here the case when at least two eigenvalues pi are different. 
Since the r.h.s. of (4.3) is the same for any i, we obtain that either all <3?j = 0, or 

^ = Pt (A + ^ i) ), (4.4) 

where A = A{r) is a function that must satisfy the equation DA + Ap = —Q^pj (since $ = Ap + ^Pj) 
and, from now on, quantities with a superscript do not depend on r. Using the intermediate substitution 
A(r) = Y(r) \\ k pk, one arrives at the solution 



•J>? f drJlir-b^+A 

3 m 



(4.5) 



4 For example in five dimensions k\ = 12 (<^2^4 + 4>3(f>4 + 4>2(j>3)- If we consider, for instance, five-dimensional RT spacetimes 
(which coincide with the Schwarzschild solution plus a possible cosmological constant), characterized by pij oc 8ij, equation (3.2) 
implies $ij oc Sij, which clearly leads to ki 7^ 0. 
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For subsequent discussions it is useful to rewrite (4.5) using partial fraction decomposition, i.e., 

u: -m,ax /-> v ' 



Pn-2 



Uk(r-b k ) Uk(r-h) 



j-y "max 

a° + ^ Fn ~\ , =a°+ £ Y, 

K=l \/Loa>K 



'(0 



(4.6) 



(K) 

where P}~ denotes a polynomial of order k in r, c)^ do not depend on r and a« denotes the multiplicity of 
bi, with a max being the maximal multiplicity. In particular, the term a can be determined by looking at the 
leading term of (4.5) in the limit r — > oo, i.e., 



so that 



°= — ^y^i 



In the case = 0, a vanishes and one has simply A = A°/Y\ k (r — &&). 

Next, the equation for Wij (see (1.12)) can be written as (but recall = 0) 



D 



D 



Wi 



P(i) ~ P(j) J P{i) ~ P(j) 



(Pj ^Pi), 



A*. 



P(i) 



(.Pj = Pi)- 



Using the above expression for $j, in both cases one can write the solution as 

1 



Wi, 



2 J drA + ($° + <S> j)r + W % 



11 ir - 6 (i) )(r - b {j) ) " 
Now, imposing — 2<3?j = ^ . (recall (3.7)) we obtain the constraint 



-2(A + $?) = ^ — 2 J drA + ($° + $°)r + Wg 



For r — >■ oo, using (4.6) and (4.8) at the leading order (4.12) implies that all coincide, i.e., 
A thus becomes 

^ = -/ + ^°n^7r> 



where we have used Ylj T\iytj( r ~ h) = P ) W\k( r ~ so that 



<F; 



r - bi r -b k 



" (r-6(i))(r-6cj)) 
The constraint (4.12) can thus be written as 



11 r - b k l \ 



(i = 2,...,n-l) 



fc n jj^i (j) L " k 

It is now useful to discuss separately various subcases with different multiplicity of eigenvalues. 



(4.7) 
(4.8) 



(4.9) 
(4.10) 

(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 



(4.16) 
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4.1.1 All the bj coincide: shearfree congruences (Robinson- TVautman spacetimes) 

In the RT case all the bi of (4.2) coincide (and can be set to zero by shifting r, if desired) and from (4.16) we 
simply obtain J2j ^ij = 0- 

On the other hand, let us assume in the following that not all bi coincide, i.e., at least two of these are 
distinct, say 62 / &3- 



4.1.2 Shearing case with all bi distinct: not permitted 

First, if all hi are distinct, i.e., a max = 1 in (4.6), we can compute explicitly the required integral using partial 
fraction decomposition 

f drTT - J— = Y 1 _ r ln(r / ~ bk \ - (b k all distinct). (4.17) 

It is thus clear that the singularity structure of the l.h.s. and the r.h.s. of (4.16) (for any i) cannot be the same 
unless ^4° = 0. Therefore, from (4.16) also X^j( r ~~ = ^- However, this condition implies that bi 

cannot be all distinct, so that this case is in fact not permitted (unless Wfj = 0, so that Wij = and therefore 
also $jj = 0, contrary to our assumptions here - however, we will see in sections 4.2 and 4.3 below that the 
case with all distinct bi is ruled out also for W-j = 0). 

4.1.3 Shearing case with at least one bi repeated 

Let us now consider the case when at least one bi is repeated, say 62, with multiplicity 1 < Q2 < n — 2, and 
there exists 63 7^ 62 (this is not a restriction since we are excluding here the RT case, which corresponds to 
02 = n — 2). 

Using partial fraction decomposition similarly as in (4.6), the left and right hand sides of (4.16) take the 
form (for i = 3) 



lhs = -2 A 



T) 

' " : + (terms with lower order poles at 62 and with no poles at 62) 



_(r-b 2 ) a2 



rhs = -2A°, , , + (terms with lower order poles at 62 and with no poles at 62), 

(a 2 - l)(r — b2) a2 

respectively, where p a2 / 0. Comparing the highest order terms in l/(r — 6 2 ) we thus get 

j4 o = Q y ^ = . (4.18) 

In particular, it follows from (4.15) that 

$i = 0, (4.19) 
so that the Weyl type is Il(abd), and (4.11) thus gives 

(r-bi)(r-bj) 

with = W$i and W$ ){i) = 0. 

Now, because of (4.20) the condition ^ ■ Wij = constraints the possible multiplicities of bi (i.e., of the 
eigenvalues of pij). First, if bi are all distinct, we immediately get Wfj = 0, as already discussed above. Similarly, 
recalling = it is easy to see that for Wfj 7^ the structures {a, a, ci, . . . , c ra _4} and {a, a, a, c\, . . . , c n -^} 



11 



are also forbidden. 5 However, the structure {a, a, b, b, ci, . . . , c n _6} is compatible with such constraints (take, 
e.g., &2 = bs, 64 = 65 and W® 4 = — W® 5 = W® 5 = —W^ / 0, and the remaining Wij equal to zero). 

To summarize, we have seen above that in type II Einstein spacetimes with a non-degenerate non-twisting 
mWAND £, one has = and 

1. if <E>jj / 0, then I must be shearfree and the corresponding spacetimes belong to the RT class. This 
includes, in particular, the result of [23] for asymptotically flat type II vacuum spacetimes (restricted to 
non-twisting case) 6 

2. if 5>jj = (type Il(abd)) the structure of eigenvalues of pij can be more generic, however for <E>jj = / Wij 
it must have the multiplicities {a, a, b, b, c\, . . . , c n -%} (or more special). 

Let us now discuss the remaining cases (in which, without loss of generality, we could assume <&ij = = Wij 
- however this will actually not be used to prove the following results). 

4.2 CaseSg^O 

Here we assume that a component of <&ijki with precisely three distinct values of i,j,k,l is non- vanishing, i.e., 
for some i ^ k, we have &i(j)k(j) / 0. This case is possible only for n > 6 because we need at least three 
possible distinct values for i, j, . . . and because of the tracefree condition <&ijik = with k / j (however, for 
n = 4, 5 necessarily ^ / and from the above discussion one already knows that RT is the only possibility). 
Eq. (2.1) gives 

$(i)U)(i)(k)(p(j) ~ P(k)) = (k^ j). (4.21) 

Without loss of generality, we can assume $2324 / 0, so that necessarily 63 = 64. 
Now, the off-diagonal (i / k) components of (1.12) read 

- D $i(j)k(j) = $i(j)k(j)P(j) + $i(j)k(j)P(k)- (4.22) 
Taking into account p^ = (r — 6fc) _1 , we arrive at 

Then, the condition Qijkj = 0, k 7^ i, implies that at least another component < I > (^)3(i)4 must be non-zero, 
say <3?5354 7^ 0, and that 65 = 6 2 . The eigenvalue structure must therefore again be {a, a, b, b, c±, . . . , c n _@} ( or 
more special). 

4.3 CaseSgUo 

Here we assume that there exists a non-vanishing component of <&ijki with all values of i, j, k, I being distinct. 
This case is possible only for n > 6. 



5 To see that {a, a, a, C2, - - - , c n -r,} cannot occur, assume 62 = &3 = &4 and all remaining eigenvalues are single. Then £^ Wij = 
with (4.20) gives W- 2 + W° 3 + Wf 4 = and = / 2,3,4). Since = this implies that the only possible non-zero 
components are W23, W24 and W34, with the conditions W23 + W24 = 0, W32 + W34 = and W42 + W43 = 0. However, this system 
admits only the solution W23 = W24 = W34 = 0, so that Wij = 0. 

6 Ref. [23] used an expansion method. Instead of the condition <f>ij 7^ 0, a condition on the asymptotic fall-off behaviour of the 
Weyl tensor was assumed there. In the present notation, that amounts to taking Wfj = = ®\fll = 0j which requires that 

$ij / (otherwise the type would be III). The assumptions of [23] in the non-twisting case were thus stronger than those used 
here (note indeed that RT spacetimes with W° / do exist [25]). (The symbols an d ®\fll are defined in sections 4.2 and 4.3 

below.) 
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(for i,j,k,l all distinct). (4.24) 
bi (or bi = bi and bj = bk) and the eigenvalue structure 

4.4 Summary 

No further cases are possible, since = ^jjli = Wij = implies the type III. To summarize, in sections 4.1- 
4.3 we have thus shown that 

Proposition 4.1. For type II Einstein spacetimes, the existence of a non-twisting, non- degenerate double 
WAND implies that the algebraic type is necessarily 11(d), i.e. <3?^ = 0, or more special, and 

1. if&ij 7^ the spacetime is shearfree ( Robins on- Trautman) , i.e. the eigenvalue structure of pij is {a, . . . ,a} 
(with a ^ 0). This is the only possibility when n = 4, 5. 7 One has &ij = A°5ij/(r — &o) n_1 > so that the 
type is D(bd) (one can set bo = by a shift of the affine parameter r). 

2. if &ij = (type Il(abd)) the structure is {a,a,b,b,c\, . . . , c n __6} (or more special; a,b,c a / 0). In 
particular, in six dimensions this means {a,a,b,b} (see also section 8). 

This is nothing but a more detailed version of points (ii) and (in the subcase det p ^ 0) (iii) of theorem 1.1 
of section 1, which are thus proven (the proof of (iii) will be completed in section 5 for any value of rank(p), 
cf. remark 5.1). In particular, (for n > 5) there are always at least two double eigenvalues. In this sense, we 
will see that the situation is different in the degenerate case. Note that spacetimes of point 1. (all explicitly 
known [25]) obey the optical constraint, whereas those of point 2. do not, in general. Examples of the latter 
in n > 7 dimensions will be provided in appendix F. 



Similarly as in the previous case, eq. (1.12) implies 

rf> — ®ijkl 
®ijki — 



( r - b (k))(r -b(i))'' 

However, since <&ijki = &kiij we obtain bi = bk and bj 
is again {a, a, b, b, c±, . . . , c h ^q} (or more special). 



5 Non-twisting I: degenerate p (0 < m < n — 2) 

Let m denote the rank of p^. The value m = n — 2 corresponds to the previously considered non-degenerate 
case, while m = defines Kundt. Therefore here we shall restrict to < m < w - 2. We need to define two 
types of indices to distinguish between non- vanishing and vanishing eigenvalues, namely 

P P = — 7 (p,q,o,t = 2,...,m + l), (5.1) 

V Op 

p z = (z,v,w,y = m + 2, ...,n-l). (5.2) 

Recall that (point (i) of theorem 1.1) in the non-twisting case with $4 ^ one has m = 0,2 and = 
diag(a, a, 0, . . . , 0) for any n > 4, so that this case does not require further investigation. In all remaining cases 
we thus have 

^ = (m/0,2). (5.3) 

In the following we will give the r-dependence of the Weyl components, which is then used to constraint the 
possible forms of p. In particular, we shall explore under what conditions all the eigenvalues of p can be distinct 
(which is not permitted in the non-degenerate case). We shall also discuss some special cases and construct 
explicit examples. 

7 Since type Il(abd) coincides with type III in those dimensions. 
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5.1 <&jj and components 

Proceeding similarly as in section 4 but (thanks to (5.2)) with the additional equation D& z = 0, one finds 



* p = Pp (A + $° {p) ), $ 2 = $°, 
1 



W pq = 



(r - b {p) ){r - b {q) ) _' 



2 J drA + ($£ + $°)r + W, 



r — b 



(p) 



where 



r -b 



E*5 / drn^-M-E^ / drII(r-6 ff ) + A 



P q^p f/ 

Now, imposing — 2<3? z = J2 P W zp + we obtain the constraints 



r-b n 



0. 



Next, imposing -2<& q = J2 P W qp + ^ 2 W ?2 gives 

-2^ + ^) = E - A- f 2 / dr ^ + + $ > + w p<\ + E^° r + ^ 



By comparing the leading terms of the l.h.s. and of the r.h.s. for r — > 00 one obtains 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



(5.8) 



(5.9) 



(5.10) 



5 - 2 $ S and components 



To complete the description of the Weyl tensor, we now give the general form of the 3>j^ and components 
(recall that these are non-zero only for n > 6). In particular, we also discuss constraints following from the 
assumption that all non- vanishing eigenvalues of pij are distinct (useful for later analysis). 

5.2.1 components 

Eq. (2.1) gives 

$ (o)(p)(o)(q)(P( ? ) - P(p)) = = ®(z)(j>)(z)(q){p(q) ~ P(p)), $(q)p(<?)z = = $ (o)p(o)z ' (p q, V ^ z). (5.11) 

Non- vanishing components must satisfy the tracefree conditions (recall that <3?jj is diagonal in the frame we are 
using) 

$poqo + %zqz = 0, ^zjwp + ^tuutu = (p / <?, V / z). (5.12) 

The r-dependence, following from (1.12), is 



(p)«(p)f 



^(p)<?(p)o 



*0 

(p)«(p)f 



r — b 



(p) 



(p)«?(p)(o) 



(r - 6(p))(r - 6 (o) ) 



®{z)v{z)w = ®°(z)v(z)w ( v ^ w )> 



r — b 



(«) 



(5.13) 
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From (5.12), the quantities in (5.13) must satisfy the constraints 



*S*. = °. Et^ = ° (p^9). ( 5 - 14 ) 



5 (o) 

*2™ = 0, E^T M = ° ( 5 - 15 ) 



(p) 



Note that, by the symmetries of the Weyl tensor, from the third equation of (5.13) it follows that &( p )q( p )o 7^ 
=>■ b q = b Q (for q / o), in agreement with (5.11); using also the second equation of (5.14) we get 

* { $ q (p)o ^°^P = dia g(«' a,b,b, Cl ...) (a, 6/0). (5.16) 

Further, the fourth equation of (5.13) gives ®( z )p(z)q / =>• b q = 6 p (for g / p), i.e., the structure is {a, a, c\ . . . } 
with a / 0, in agreement with (5.11). For <3? 2 ( p )t,( p ) / 0, the second equation of (5.15) also implies {a, a, c\ . . . } 
with a^0. 

We are interested, in particular, in determining what are the necessary conditions in order to have all the 
non-vanishing eigenvalues distinct. From the above observations it follows that one necessarily has $(p) g (p) = 

®(z)p(z)q = ®z(p)v(p) = 0) therefore the only non-zero components can be the $ z ( w ) v ( w ), however with the 
constraint Q Q ZWVW = (eq. (5.15)). It is easy to see that for this to be satisfied in a non-trivial way, indices 
z, v . . . must run at least over four values, i.e. there must be at least four zero eigenvalues of pij (excluding 
Kundt, the spacetime must thus be at least seven-dimensional). In other words: 

• all non-zero pij are distinct and m > n — 6 =>• = 0. 
5.2.2 components 

From (1.12) and the symmetries of the Weyl tensor one finds 

pq(o)(t) _ pz(q) W o fC - ir) 

(r-6( ))(r-6( t )) r - 6 (?) y 

^pqoz 0, ^pgzto 0, ^jjotuj 0- (5.18) 

Similarly as in the case of 3>^, we observe that if ^ pqo t / then 6 t = 6 g and b Q = b p (or b t = b p and 6 Q = b q ), 
i.e., 

$2 t /0^p = diag(a,a,6,6, Cl ...) (a, 6/0). (5.19) 

Additionally, if & pzqw / then b q = 6 P , so that the structure is {a, a, c\ . . . } with a / 0. 

Having all the non-vanishing eigenvalues distinct thus requires Q pqo t = = <& pzqw and the only non-zero 
components can be $ zvwy . Therefore we conclude again that indices z,v ... must run at least over four values 
unless ^jjli = (and, again, excluding Kundt, the spacetime must thus be at least seven-dimensional), i.e., 

• all non-zero pij are distinct and m > n — 6 =>■ ^\fli = 0- 

5.3 Case $^ = (type Il(abd)) 

This case (non-trivial only for n > 5) is obtained by setting 

A = = K°, $° = (5.20) 
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in the results obtained in 5.1. For the W^ components one thus has 

W m = l 1 v I V W p* = — T^' W gv = W? v , (5.21) 

(r - b (p) )(r - b (q) ) " r-b {p) 

with 

W° W° 

EK = °, E^ = ' = ?^t = a <5 - 22) 



In view of these constraints (and of the properties of Wij) we can briefly comment on some special cases: 

• <Z>ij = 0, m = 1 =► W pq = = W pz , 

• = 0, m = 2 or m = 3 => W^g = 0, 

• = 0, m = n - 3 =>■ W z „ = = VF p2 , 

• 3>jj = 0, m = n — 4 or m = n — 5 =>■ W a „ = 0. 

For special values of n and m some of these can hold simultaneously, thus leading to Wij = 0. Recalling 
the trivial implications m = =>■ Wp ? = = PF P 2 and m = n — 2 =>■ W^t, = = Wp Z (valid also for $y ^ 0), 
we have in particular 

• for n = 4 and m = 0, 1, 2, = = 0, 

• for n = 5 and m = 0, 1, 2, 3, $y = = 0, 

• for n = 6 and m = 1, 3, <3?jj = =>• Wjj = 0. 

While the first two implications simply reproduce the known result that $ij = <^ = for n = 4, 5 (for 
any permitted m), the last remark will be useful for later purposes. 

Using (5.21), (5.22) and a reasoning similar to that of section 4.1.3 (in the paragraph after eq. (4.20), 
including footnote 5) one can also show that 

Wp q ^0^ p = di&g{a,a,b,b,c 1 ...) (a, 6/0), (5.23) 

which will be useful in the following. 

5.3.1 No repeated non-zero eigenvalues 

In case all the b p are distinct, equations (5.22) imply 

W% = = W° q , (5.24) 

so that the only non-zero components of Wij are the W zv = W zv , with Ylv^zv = ®- This implies that, in 
order to have Wij / 0, one needs that indices z,v,. . . run at least over four values, i.e., the cases m = n — 3, 
m = n — 4 and m = n — 5 imply Wij = 0. Thus we have proven: <I>jj = 0, all non-zero pij are distinct and 
m > n — 6 =>• Wij = 0. Together with the results of subsection 5.2, we see that if such assumptions hold all 
boost weight zero components must actually vanish, so that we can conclude 

Proposition 5.1. For type Il(abd) ($ij = 0) non-twisting Einstein spacetimes with (degenerate) p of rank 
n — 6 < m < n — 2, the eigenvalue structure of p is {a, a, c± . . . , 0, . . .} (with a / 0) or more special (i.e., at 
least two non-zero eigenvalues of p coincide). 

The condition on m means that either m = n — 3 or m = n — 4 or m = n — 5. We will see in section 5.4 below 
that in the case m = n — 3 the assumption <3?jj = can in fact be dropped (since $jj 7^ is not possible if all 
non-zero eigenvalues are distinct). For m = n — 4 the eigenvalue structure is thus {a, a, 0, 0, c\, . . . , c n -§} and for 
m = n — 5 it is {a, a, 0, 0, 0, C\ , . . . , c n -j} (while for any smaller m there are, of course, n — 2 — m > 4 vanishing 
eigenvalues). Recall also that for m = n — 2 (non-degenerate case) we had a similar result, cf. proposition 4.1. 
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5.4 The special case m = n — 3 

For n = 4 this case is not possible since it gives m = 1 (see section 5.5). In higher dimensions, the simplest 
spacetimes of this class are obtained by taking a Brinkmann warp [29] (see also [30]) of a type II Einstein 
spacetime possessing a non-degenerate mWAND with a single (spacelike) extra-dimension (cf . also [30] ) . They 
include, e.g., static black strings and this case may thus be of special interest. 

Since there is only one vanishing eigenvalue of pij, i.e., z can take only one value, we have W zv = and, 
from (5.8), 

^ W° 



r -b„ 
p p 



Therefore (5.10) reduces to 



(n-3)^ + ^ = ^$°. (5.26) 



{3} 

From the results of subsection 5.2, when m = n — 3 the only possible non-zero components of can be 
®{p)q{ p )oi with q / o (in particular, & pzqz = ® p ( z )q(z) = 0> see the first of (5.14)), which satisfy 

E%^ = ° (p*9)> ( 5 - 27 ) 



o 



r - 0( o ) 



while the only possible non-zero components of are ^pqot (^pgoz = 0, see the first of (5.18)) . From 

sections 5.2.1 and 5.2.2 (see (5.16), (5.19)) it thus follows that if fcjjj, / («■ $(pW p)o / here) or / 
(•w> ^p^of 7^ here) then the eigenvalue structure is {a, a, 6, b, c\ . . . }, where a, 6 7^ 0. Moreover, it is easy to 
see that when m = n — 3 both and can be non-zero only for n > 7. This is of interest only in six 

dimensions (^^ = = ^jfli identically for n = 4,5) and can be summarized as 
• n = 6,m = 3=>*W =0 = *£>. 

If we additionally assume <J>ij = (then also = = Wp 2 , see comments after (5.22)), then some of the 
components W pq , ^^l^^ or <&*!pq 0t have to be non- vanishing (otherwise the spacetime would be of type III). 
Then using also results of section 5.3 (eq. (5.23)) it follows that type II spacetimes with m = n — 3 and <frjj = 
must have the eigenvalue structure {a, a, b,b,c\ . . . c ra _7, 0}. 

5.4.1 No repeated non-zero eigenvalues 

If all bp are distinct from the second of (5.25) we obtain 

W% = 0, (5.28) 

so that W pz = and all the <3?° coincide, namely, 
Therefore 

A=-f° + A°T\ _J_ , (5.30) 
±± r — b v 

p p 

and one can basically proceed as in the non-degenerate case. Comparing the structures of the poles of the l.h.s. 
and r.h.s. of (5.9) one again arrives at 

,4° = 0, W® q = 0. (5.31) 
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Combining this with the previous results we conclude & g = 0, <E> 2 = 0, W qz = 0, W gp = 0, W zv = 0, i.e. 

$ ij = = W ij . (5.32) 

Since here m = n — 3 > n — 6, we can use the results of subsection 5.2 to conclude that also = = 3>|^, 
so that all b.w. components of the Weyl tensor vanish. Thus we have proven: all non-zero p^ are distinct 
and m = n — 3 =4> all b.w. components vanish. In other words, in the case m = n — 3, the eigenvalue structure 
of p must be {a, cl, c\, . . . , c n _5, 0} (or more special) . If, additionally, *&ij — 0, there is in fact another repeated 
non-zero eigenvalue [20] 8 (see the comment just before section 5.4.1), so that Proposition 5.1 can accordingly 
be reformulated as (see also [20]) 

Proposition 5.2. In a type II Einstein spacetime of dimension n > 4 with a non-twisting multiple WAND of 
rank m = n — 3, the eigenvalue structure of p is {a, a,c\,. . . c ra _5, 0} where a, c a / (and necessarily ^f- = 
if n > 5). If additionally, <3?jj = 0, the structure becomes {a, a,b,b,c±, . . . c n -j, 0} where a, b, c a / 0. 

Remark 5.1. Together with propositions 4.1 and 5.1, this shows that, for any possible value of m (0 < m < n— 2), 
the mWAND of type Il(abd) (i.e., <&ij = 0) Einstein spacetimes has p of the form {a, a, b, b, c\, . . . c n -o} (where 
a/0ifm = n- 5orm = n- 4 and a, 6^0 if m = n — 3 or m = n — 2). In particular, point (hi) of 
theorem 1.1 is thus now also proven. 

Remark 5.2. In five dimensions (for which always <&ij ^ 0), recalling also the results for the non-degenerate case 
(proposition 4.1), this enables us to list all the admitted eigenvalues structures, for m = 3, 2, 1, 0, respectively: 
Pij = diag(a, a, a), = diag(a, a, 0), p^ = diag(a,0, 0), or pij = (Kundt). See [12] for more details and 
examples. 

Remark 5.3. In six dimensions, we conclude that the case m = 3 requires p^ = diag(a, a, b, 0). However, as 
shown above in this case one has = = and additionally the implication <3?jj = =^ Wij = holds 
(see subsection 5.3). Therefore one necessarily has <J>f,- / when n = 6 and m = 3 (or all Weyl boost zero 
components would vanish). See section 8 for more details. 



5.5 The special case m = 1 

From a viewpoint complementary to that of the previous subsection, it is also interesting to analyze the case 
when only one eigenvalue is non-zero. Let us associate the index p = 2 to this eigenvalue, while z, v = 
3, ... ,n — 1. Using the general equations presented above, one obtains (cf. also [28]) 

<I>, = $o, <D 2 = -I^d>o, W 2z = <S>° z , £< = X (5.33) 

Z V 

while clearly here W pq = (note that for n = 4 one gets $° = and $ij = = Wij, so this case does not 
occur). Therefore all components of and Wij are r-independent. Note that <3?jj = <^ = 0. When this 
happens the only non-zero components of Wij can be the W® v , but with the constraint J2 V ^zv = 0- Therefore 
it is easy to see that Wij can be non-zero only if z, v can take at least four values, i.e., z, v = 3,4,5,6, . . ., so 
that n > 7. To summarize 

• n < 7, m = 1, ®ij = W^ = . 

This can be understood as a special instance of the results of subsection 5.3. This result is of interest in six 
dimensions, since one always has $ij = =^ Wy = in four and five dimensions. Namely, for n = 6 the case 
{a, 0, 0, 0} is forbidden if § tj = / W t j. 

8 We thank Lode Wylleman for pointing this out. 
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From the results of subsection 5.2, when m = 1 the only possible non-zero components of can be 
^(z)v(z)w (with w/to), which satisfy 

C^ = (v^w), (5.34) 

while the only possible non-zero components of are & Z vwy It is easy to see that, since m = 1, both sets 
of components vanish identically unless n > 7, so that 

• n = 6, ro = 1 =► *g> = = *£>. 

Remark 5.4. Since for n = 6, m = 1 one has = =^> Wjj = (see above), we conclude that in six dimensions 
spacetimes with m = 1 require <&jj / 0. 

5.5.1 Examples 

For any n > 5, simple examples of spacetimes with m = 1 (and thus $4 = 0) are given by: 

1. dS 3 xS n _ 3 (or AdS 3 xH n _ 3 ): n > 5, ^ + 0, A ^ 0, 

2. as above but with a Brinkmann warp: n > 6, / 0, also A = is possible, 

3. Minkowsky x(Ricci-flat) n _ 3 : n > 7, $y = 0, A = 0. 

In the first case, if we restrict to n > 7 the second factor space can be in fact any Einstein space (with Ricci 
scalar given by 2(n — S)K). The metric of the first and third case can thus be given in a unified way as 

K \~ 2 



ds 2 =[l + —urj (2d«dr + r 2 dx 2 ) + dS 2 , (5.35) 

where dE 2 is the metric of an (n — 3)-dimensional Einstein space with Ricci scalar = 2(n — 3)K and K is a 
constant related to the n-dimensional Ricci scalar by R = 2nK. A geodesic twistfree mWAND £ (degenerate 
but expanding along x) is 

£ a dx a = du. (5.36) 

Note that here r is an affine parameter along t only in the case K = 0. Such product spaces are of type D [15] 
and in fact any null vector field tangent to the three-dimensional Lorentzian factor is an mWAND [15, 16]. 
Using these results metrics for the second case can also be constructed straightforwardly. These will be also of 
type D [30]. 

5.6 Case when all non-zero eigenvalues are distinct 

By combining the results of subsections 5.3 and 5.2 we have 

• all non-zero pij are distinct, m = n — 4 or m = n — 5 and $ij = all b.w. components vanish, 

• all non-zero pij are distinct and m = n — 3 all b.w. components vanish. 9 

Therefore type II Einstein spacetimes for which all non-zero pij are distinct require that at least one of the 
following holds (obviously m > or we would simply have Kundt): 

1. 1 < m < n — 6 (i.e., at least four eigenvalues of pij vanish), 

2. $jj / and 1 < m < n — 4 (i.e., at least two eigenvalues of p^ vanish). 



3 By adding the results of [8] we can in fact conclude that in both cases the whole Weyl tensor must vanish. 



19 



Case 


n 


Possible m 


Examples 


*ij + o 


5 


1 






\) 


1 2 






7 








8 


1,2,3,4 






n 


1,2,3. . . , n - 4 


( r\ 1 (] l> «i i oXh« O I H ,, o ™ 1 

V J V w<~r^ z— m V n — z— m J 


*ij = o 


7 


1 






8 


1,2 






n 


1,2,3,. . . , n — 6 


Mink m+2 x (Ricci-flat) n _ 2 - m 



Table 1: Non- twisting Einstein spacetimes for which there are no repeated non-zero eigenvalues of pif permitted 
values of the rank m of pij in various dimensions (necessarily 1 < m < n — 4 thanks to propositions 4.1 and 5.2). 
As discussed in the text, different examples can also be generated using Brinkmann's warp, or by replacing the 
S n _2- m (H n _2-m) factor by a more general Einstein space. Note that here solutions with m > 1 violate the 
optical constraint. 



The first case can occur only for n > 7: for n = 7 one can have only pij = diag(a, 0, 0, 0, 0), for n = 8 either 
Pij = diag(a, 0, 0, 0, 0, 0) or pij = diag(a, b, 0, 0, 0, 0), etc.. The second case can occur for n > 5: for n = 5 one 
has simply pij = diag(a,0, 0), for n = 6 either pij = diag(a, 0, 0, 0) or pij = diag(a, b, 0, 0), for n = 7 either 
Pij = diag(a, 0, 0, 0, 0) or pij = diag(a, b, 0, 0, 0) or pij = diag(a, b, c, 0, 0, 0), etc.. These results are summarized 
in table 1. For such solutions the optical constraint is clearly violated if m > 1. The previously discussed case 
m = 1 has been included for completeness but it is of course "trivial" in this context since there is a single 
non-zero eigenvalue. 

5.6.1 Examples 

The case n = 4 is not possible. Examples with m = 1 have been already presented in subsection 5.5 and will 
not be discussed again here. For m > 2 the n = 5 case is also forbidden (cf. Remark 5.2 and [12]). For n = 6 
the only possibility is m = 2 (with $^ ^ 0) and one can construct explicit Einstein spacetimes by taking 
dS4xS 2 or AdS4xH 2 . These can then be extended to any dimension n > 7 by a simple Brinkmann warp (and 
if one starts from dS4xS 2 the cosmological constant of the resulting spacetime can be arbitrary, cf. [30]), so to 
have explicit examples of the type pij = diag(a, b, 0, . . . , 0) (i.e., with m = 2) for any n > 5. 

However, according to the previous comments (see also table 1), for n = 7 one can have m = 2,3, for 
n = 8 it is possible m = 2,3,4 and so on: in general, and for any n>5, 2 < m < n — 4 explicit metrics can 
be constructed similarly as in six dimensions by taking dS m +2 x S n _2-m or AdS m +2 x H Tt _2-m (where the two 
factor spaces must have Ricci scalars given by (m + 2)(m + 1)K and (n — 2 — m)(m + 1)K, respectively), and 
appropriate hypersurface orthogonal null congruences living in the (A)dS m+ 2 factor (explicit examples are given 
below). These have necessarily <&jj / and A / 0. Using Brinkmann's warp one can also generate Ricci- flat 
solutions (again with ^ 0) for which, however, the stronger restriction n>6, 2 < m < n — 4 holds (at least 
if this method is used). Products of Ricci-flat spaces are possible if each factor is at least four-dimensional 
(therefore there are always at least four vanishing eigenvalues) and the Lorentzian factor is flat and fall into 
the &ij = class, with A = (these metrics belong to the type D pp -waves mentioned in [21]). 

Summary: 

1. dS m+2 x (or AdS m+2 xH n _2-m) 10 : n > 5, 2 < m < n - 4, $^ / 0, A / 0, 

10 Note that if n — 2 — m > 4 the second factor space can be in fact any Einstein space with Ricci scalar given by (n— 2— m)(m+l)K. 
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2. as above + Brinkmann, also A = 0: n > 6, 2 < m < n — 5, <J>jj / 0, 

3. Minkowski m+ 2x(Ricci-nat) n _ 2 - m : n > 7, 2 < m < n — 6, <&jj = 0, A = 0. 

Note that all the above spacetimes also belong to the Kundt class, although w.r.t. a null congruence different 
from the one which is of interest for our discussion (which is expanding and shearing). In fact they all admit 
oo m mWANDs: similarly as in section 5.5.1, they are of type D (with = 0) and any null vector field tangent 
to the Lorentzian factor is an mWAND - so one can always also find an mWAND with an optical matrix 
having all non-zero eigenvalues identical. However, this mWAND "degeneracy" is not the generic situation 
(see appendix F for an example of a class of type D spacetimes which admit (only) two double WANDs, both 
violating the optical constraint). 

Explicitly, for cases 1 and 3 one can take, e.g., the metric 

ds 2 = Q 2 (2dudr + Pgdx 2 q ) + d£ 2 , (5.37) 

where q = 2, . . . ,m + 1, d£ 2 is the metric of an (n — 2 — m)-dimensional Einstein space with Ricci scalar 
i? s = ( n - 2 - m){m + l)K and 



P q =r-b q , fT 1 = 1 + — 



K 
I 



(5.38) 



in which b q and K are constants. The first factor space in (5.37) is a space of constant curvature with Ricci 
scalar given by (m + 2)(m + l)K, so that the Ricci scalar of the full spacetime is 

R = n(m + \)K. (5.39) 

The geodesic, twistfree mWAND £ and the remaining frame vectors are given by 

£ a dx a = du, n a dx a = Q 2 dr, m^ q)a dx a = UP^dx^, (5.40) 

while vectors m^ a will depend on the specific form of dS 2 . The optical matrix pij is given by 

Ppq = 6 (p)q 3 P ' ' Ppv = = p vz . (5.41) 

In the case K = (i.e., A = 0) one thus has p pq = S^ q /(r — fo( p )). Note, however, that r is not an affine 
parameter when K ^ (one has £ = Q,~ 2 d r ). 

For case 2 one can use the method illustrated in detail in [30]. 



6 Counterexamples 

In [12], a counterexample to the converse of the five-dimensional "shear-free" part of the GS theorem was 
presented, thus demonstrating explicitly that the condition that the optical matrix p admits a canonical form 
compatible with a geodetic mWAND is not sufficient for the null geodetic being an mWAND. Similarly, also 
our results above give conditions that are (necessary but) not sufficient. In order to demonstrate that, here we 
present a few "counterexamples", i.e., certain Einstein spacetimes that admit a non- twisting null (thus geodesic) 
vector field £ with p taking one of the permitted "canonical forms" (cf. theorem 1.1 and propositions 4.1, 5.1 
and 5.2) and yet with £ not being an mWAND. Note that such counterexamples will necessarily be shearing: 
a twistfree shearfree null vector field is automatically an mWAND [10,25]. 
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6.1 Non-degenerate p = {a, a, b, b, ci, . . . , c„_ 6 } 



In the non-degenerate case a "counterexample" to the canonical form {a, a, b, b, ci, . . . , c n -§} (with a, b,c a 7^ 
and 0/6 - see points (ii) and (iii) of theorem 1.1 and proposition 4.1) is given by the following six-dimensional 
Ricci-flat spacetime, which belongs to a class of metrics considered by Robinson (as described in [31]), 

ds 2 = 2dudr + 2 cosh 2 rdwdw + 2 sin 2 rd(d(. (6.1) 

Here the hypersuface orthogonal null vector field £ a dx a = du is not an mWAND (not even a single one) and 
yet the corresponding p has the eigenvalue structure {a, a, b, b} (with a 7^ 7^ b, a 7^ b). Associated to I there 
is also an optical structure [31]. Note, however, that spacetime (6.1) is a type N pp-wave (albeit in so-called 
Rosen coordinates, see, e.g., metric (39) of [29]), with a covariantly constant mWAND given by n a dx a = dr. 



6.2 Case n — 6<m<n — 2 

Direct products of (6.1) with flat extra dimensions clearly provide non-WANDs with eigenvalues {a, a, b, b, 0, . . . , 0}. 
In particular, for n = 7 one obtains {a, a, b, 6,0}, which is thus a counterexample to the p-form of Proposi- 
tion 5.2. For n = 8, 9 this is also a counterexample to the p-form of Proposition 5.1. Another counterexample 
to Proposition 5.2 can be obtained for n = 6 by taking the direct product of a vacuum black ring [32] with a 
flat dimension: in the coordinates of [33], the hypersurface orthogonal null vector field (also considered in 5d 
in [19] for different purposes) 

Ldx a = ^^dy + dV, (6.2) 

corresponds to p of the form {a, a, b, 0} and is not a WAND. Alternatively, the same form {a, a, b, 0} can be 
obtained by taking the direct product of a 5d vacuum static KK bubble [16] with a flat dimension (with a 
non-WAND £ a dx a = dt + V(r)~ l l 2 dr - this 6d spacetime is of type D but the mWANDs are different from 
I [30]). 



6.3 Case when all non-zero eigenvalues are distinct 

Examples of spacetimes with a non-twisting, non-WAND I having all non-zero eigenvalues distinct can be 
constructed as a direct product of a 4d pp-wave with flat space, i.e., 

ds 2 = 2dudr + cosh 2 rdx 2 + sin 2 rdy 2 + dzidz 1 . (6.3) 

Here, £ a dx a = du has a p with eigenvalue structure {a, b, 0, . . . , 0} (with a 7^ 7^ b, a 7^ b). This is the form of 
p discussed in section 5.6 (see also table 2 for the n = 6 case). Similarly as the spacetime (6.1), also (6.3) is a 
type N pp-wave with mWAND n a dx a = dr. 



6.4 Case m = 1 

Examples of non-twisting null congruences with m = 1 (the case discussed in section 5.5) that are not multiple 
WANDs can be obtained from the (Ricci-flat) Newman- Tamburino lift considered in [12], i.e., 

ds 2 = r 2 dx 2 + x 2 dp 2 dudx - 2dndr + x~ 2 [c + ln(r 2 x 4 )] du 2 + dzidz 1 , (6.4) 

where c is a constant. In this case £ a dx a = du is a single WAND with optical matrix of the form {a, 0, ... , 0} 
(with a / 0). 
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7 Existence of totally geodesic null two-surfaces 



Theorem 1.1 comprises the main results proven in the previous sections. In particular, it implies that there 
always exists at least a pair of repeated (possibly vanishing) eigenvalues and that (recall that n > 6) there are in 
fact at least two such pairs when, e.g., det p/0 (see also proposition 4.1) or <E>jj = (see also propositions 5.1 
and 5.2). The following proposition (see also the definitions in appendix C, in particular eq. (Cll)) elucidates 
the geometrical meaning of this fact (and provides a connection with the standard 4d version of the Goldberg- 
Sachs theorem, cf. appendix C.l). 

Proposition 7.1. In an algebraically special Einstein spacetime of dimension n > 6 with a non-twisting 
mWAND £, using a parallelly transported eigenframe of p 

1. if the eigenvalue structure of p is {a, a,c±, . . . , c n _5, 0} with c a / a / then the totally null distribution 
V23 is integrable with totally geodesic integral surfaces, 

2. if the eigenvalue structure of p is {0,0, ci, ... , c„_4} with c a / then the totally null distribution V23 is 
integrable with totally geodesic integral surfaces, 

3. if the eigenvalue structure of p is {a, a,b,b,c±, . . . , c„_6} with b 7^ a and a / c a 7^ b then the totally null 
distributions V23 and D45 are integrable with totally geodesic integral surfaces (either a or b can be zero). 

Proof. The proof is similar in all cases and relies on the use of the Ricci identity (Ilk, [10]) (equivalent to the 
NP equation (A4, [11])), which for algebraically special Einstein spacetimes reduces to 

1 1 

%|ft|fe] = L i[j\Pi\k] + T ip[jk] + PuM[jk] + Pi[j\Mi\k], (7.1) 

where 5{ = m^V a (recall also the definitions in (1.6)). Since here £ is twistfree we have p^} = in the above 

equation. Assume the form of p specified in 1., i.e., P22 = P33 = a / (and all remaining eigenvalues different 

23 

from a) and consider (7.1) for i,j, k = 23k/32k (with k ^ 2,3). This gives M k3 = = M k2 - Then, (7.1) with 

2 3 

1, j, k = 22k /2,3k gives M k2 = M k3 , which completes the proof of point 1. (see eqs. (C15)). The proof of point 

2. is identical. The proof of point 3. is analogous, after extending the reasoning to the index pair (45); in this 
case, however, either a or b can be zero, provided all c a / 0. □ 

Remark 7.1. In the case there are further pairs of repeated eigenvalues, the corresponding two-spaces will also 
define integrable totally geodesic distributions, provided these pairs are not repeated. In particular in even 
dimensions, if there are (n — 2)/2 (non-repeated) pairs of equal eigenvalues there will be (n — 2)/2 such totally 
geodesic two-surfaces. 

Remark 7.2. Point 1. of the above proposition includes, in particular, all non-twisting but expanding type 
III/N Einstein spacetimes (for which c a = [8]) and type II non- twisting non-Kundt spacetimes with <tf- 7= 
(point (i) of theorem 1.1). Point 3. is relevant to the case detp 7= with 5>jj = (proposition 4.1 and point 
(iii) of theorem 1.1), which can only be of type Il(abd). 

Remark 7.3. The proposition does not explicitly include the n = 5 case. However, from Proposition 4 of [12] 
one can easily obtain the corresponding results (also in the presence of twist). In particular, in the non-twisting 
case T>23 is totally geodesic when the eigenvalue structure of p is either {a, a, 0} or {0, 0, a} (with a 7= in both 
cases). The same is true in the case {a, a, a} 7= 0, although this does not follow from (Ilk, [10]) (instead, one 
can easily adapt the argument of footnote 10 of [12]). 

Remark 7.4. Recall that in 4d the Goldberg-Sachs theorem implies that £ is geodesic and the eigenvalue 
structure of is {a, a}, which is equivalent to V23 being integrable (and thus automatically totally geodesic, 
in 4d, see appendix C.l). 
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8 Type II spacetimes with a non-twisting multiple WAND in six dimen- 
sions 



As already mentioned, five-dimensional algebraically special Einstein spacetimes have been studied in detail 
in [12] and, for n = 5, the results obtained above are also contained in that reference. The next lower dimension 
to consider is thus n = 6 and in this section our general results are specialized to this case. This is of interest 
also because of a qualitative difference between n < 6 and n > 6 dimensions: for n < 6 type II spacetime 
necessarily have $jj / (since type II coincides with type 11(c) for n = 4,5), while this is not so in higher 
dimensions. 

8.1 Permitted forms of p 

One has to consider the various cases m = 4, 3, 2, 1 (m = is Kundt). 

First, recall that for m / 2 point (i) of theorem 1.1 implies that <E>4 = 0. Let us thus first discuss this case. 

Now, for m = 4 proposition 4.1 tells us that if (fr^ / Owe are in the RT class, i.e. p = diag(a, a, a, a) / 
(in which case <&ij oc 5ij), otherwise we have p — diag(a, cl, b, 6) (with a, b 7^ 0). 

For m = 3 we learn from Proposition 5.2 that p = diag(a, a, b, 0), while from section 5.4 we know that 
$W = = <S>\f kl and ^ / (Remark 5.3). Substituting the form of p into (3.2)-(3.4) and (A5) further 
reveals that &ij = diag(a, —a, 0, 0) / 0. 

For m = 1 one obviously has p = diag(a, 0, 0, 0) / 0. From section 5.5 we know that = = <E>^ and 

/ (Remark 5.4). 

Finally, for m = 2 we can have two possibilities. If 7^ point (ii) of theorem 1.1 shows that p = 
diag(a, a, 0, 0) 7^ 0. On the other hand, if = = we still get p = diag(a, a, 0, 0) / from Proposition 5.1, 
whereas if $4 = / the more general form p = diag(a, b, 0, 0) is permitted (with a 7^ ^ 6). 

We have checked in the various permitted cases that the constraints (3.2)— (3.4) and (A1)-(A6) allow for 
nontrivial type II Weyl tensor. These results are summarized in table 2, together with a few examples. Examples 
of spacetimes with distinct non-zero eigenvalues follow from the general discussion of section 5.6.1, which need 
not be repeated here. 

8.2 Additional examples 

In addition to the examples mentioned in table 2 and in section 5.6.1, here we give two more solutions with 
m = 1 in six dimensions (and, necessarily, 7^ 0). 

8.2.1 Example with A / 

We can use, e.g., AdS4xH 2 in the form 

ds 2 = ^-(-dt 2 + dp 2 + p 2 d<t? + dz 2 ) + ^{dx 2 + dy 2 ). (8.1) 

For the null congruence 

£ a dx a = dt + dp, (8.2) 

with the spacelike frame vectors mp) = z/9 _1 /(v / 3/3)3 ( ^ ) , m( 3 ) = z/(y/3j3)d z , = xfi~ x d x , m( 5 ) = x^~ x d y , 
one finds pij = diag(a, 0, 0, 0), where a = z 2 /(3/3 2 p). 
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Case 


m 


Possible form of p 


o.c. 


o.s. 


I.T.G. 


Examples 


®ij + o 


4 




diag(a, a, a, a) 


V 






RT 




3 




diag(a, a, b, 0) 


if b = a 


X 


P 23 (if b / a) 


?(t) 




2 




diag(a, b, 0, 0) 


if b = a 




^45 (^23 if b = a) 


dS 4 xS 2 




2 


(n * o) 


diag(a, a, 0, 0) 


V 


V 




? 




1 


diag(a, 0,0,0) 


V 


X 




dS 3 xS 3 









diag(0, 0,0,0) 


V 






Kvmdt 


*ij = o 


4 




diag(a, a, b, b) 


if b = a 




^23, ^45 (if & + a) 


?(t) 




3 




X 








X 




2 




diag(a, a, 0, 0) 


V 




^23) ^45 


? 




1 




X 








X 









diag(0, 0,0,0) 


V 






Kvmdt 



Table 2: Permitted forms of the optical matrix associated with a non- twisting mWAND in a six-dimensional 
type II spacetime in the two cases 3>jj / and 3>jj = (subtype Il(abd)). Recall that <J>4 ^ is possible 
only for m = 2. A question mark indicates subcases that are in principle permitted but for which no examples 
are known (so we do not claim that these classes are non-empty), while X corresponds to forbidden subcases. 
Although we are not aware of examples with m = 3 and <E>jj / 0(= <E>^) having the most general structure, 
special solutions with a = b (e.g., a static black string) can be obtained by warping once a type II 5d Einstein 
metric. Similarly, generic m = 4 solutions with $jj = are not known, but in the special case a = b they 
exist within the RT family [25] (the "/i = 0" subcase). These special examples are denoted by a f in the 
table. In the O.C. column we have indicated if/under what conditions p satisfies the optical constraint. In the 
O.S. column whether V = Span{m( 2 ) + im^,m^ + imt 5 \,£} defines an optical structure and in the I.T.G. 
column we listed integrable complex two-dimensional distribution that admit totally geodesic integral surfaces 
(cf. section 8.3 and appendix C). The symbol * means that the corresponding statement is true in a special 
subcase specified in the text (see section 8.3). 



8.2.2 Example with A arbitrary 

One can alternatively construct an example with pij = diag(a, 0, 0, 0) by starting from a null geodesic congruence 
in dSsxS 2 and adding an extra dimension by Brinkmann's warp. For example, 



ds 2 = f(z) 
with the null congruence 



-V(r)dt 2 + p^ydr 2 + r 2 d* 2 + (3 2 (d9 2 + sin 2 6d<f> 2 ) 



2 



£ a dx a = dt + V^dr, (8.4) 
where 

V(r) = l-^, f(z) = -Xz 2 + 2dz + b, d 2 = ^L-\b. (8.5) 

Using the spacelike frame vectors m^) = f~ 1 ^ r ~ 1 d x , m(3) = f~ 1 ^ 2 fi~ 1 de, = f^ 1 ^ 2 ^ 1 sin -1 Od^, m( 5 ) = 
f 1 / 2 d z , one finds pij = diag(a, 0, 0, 0), where a = jj. 

It follows from [15,30] that $>jj / 0, as it should be. In general A can have any sign but, in particular, we 
can set A = by choosing the parameter A = in f(z). 

8.3 Optical constraint, optical structures and totally geodesic complex two-dimensional 
null surfaces 

First, in the non-twisting case the optical constraint (2.7) is satisfied iff the optical matrix p possesses only one 
(possibly repeated) non-zero eigenvalue plus, possibly, some vanishing eigenvalues. This explains column O.C. 
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("optical constraint") of table 2. 

Next, when I is non-twisting one can always take an eigenframe of p to be parallelly transported [24, 

i 

28] and thus in particular set Mjo = 0. It is then obvious that to each pair of repeated eigenvalues of p 
there corresponds an integrable two-dimensional totally null distribution, defined by two unit vectors of the 
corresponding eigenspace and by £ (cf. eq. (C13)). 

Further conditions on the Ricci rotation coefficients can be obtained by considering the Ricci identity (Ilk, 
[10]), i.e., eq. (7.1), similarly as in section 7. Consequences of this equation enable one to prove the statements in 
columns O.S. ("optical structure") and I.T.G. ("integrable totally geodesic" complex two-dimensional subspace) 
of table 2. Since this is rather straightforward, let us just exemplify it to show that when pij = diag(a, a, b, b) 
(with b / a) the maximally totally null distribution V = Span{m( 2 ) +im^,m^ +im^,£} defines an optical 

4 4 

structure (cf. eqs. (CIO)). Namely, from (7.1) with k = 224 and k = 334 one gets M22 = M33, while 

5 5 

k = 225 and i,j,k = 335 give M22 = M33. Similarly, by taking k = 442/552 and i,j,k = 443/553 

2 2 3 3 

one obtains, respectively, M44 = M55 and M44 = M55. Next, with i,j,k = 324,234,325,235 one finds 

3 232 5454 

M42 = M43 = M52 = M53 = 0, while i,j,k = 524,425,534,435 give M24 = M25 = M34 = M35 = 0. This 
implies that eqs. (CIO) are satisfied, as we wanted to prove. 

Let us finally remark that in the Robinson- Trautman case, corresponding to the first row of table 2, eq. (7.1) 

i 

reduces to Lu = (cf. also the results of [25,28]) and thus does not constraint the Mjk- However, in the 
special subcase (cf. [25]) in which the transverse space is of constant curvature one can show by inspection 
that there exists a natural parallelly transported frame such that D defines an optical structure and T>23, ^45 
correspond to totally geodesic null two-surfaces. 11 In fact the real mu\ can be paired arbitrarily in this case, 
giving rise also to other optical structures and totally geodesic null two-surfaces, and since there is a continuous 
freedom in spins there exist actually infinitely many optical structures. This partly answers a question raised 
in remark 5.3 of [27], by showing that also in even dimensions a spacetime may admit more than 2[™/ 2 l optical 
structures (noting also that in this example replacing I ■<->■ n gives other optical structures). 

8.4 Type III/N spacetimes 

By the results of [8] , in this case there are only two admitted structures of p: either p = (Kundt spacetimes) or 
p = diag(a, a, 0, 0) / 0. Since eq. (7.1) does not contain any Weyl components of zero boost weight, the above 
discussion still applies and implies that in the case p = diag(a, a, 0,0) / both V23 and P45 correspond to 
totally geodesic null two-surfaces and that V defines an optical structure - explicit examples can be constructed 
as explained in [24]. 
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A Algebraic constraints for type II spacetimes 
A.l General mWAND 

For completeness and future reference, let us give here some additional algebraic constraints following from the 
Bianchi and Ricci identities for type II Einstein spacetimes in arbitrary dimensions. 

11 This is a straightforward extension of the discussion in the 5d case, see footnote 10 of [12]. 
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Note that rhs of eq. (1.12) is not symmetric under ij «-> kl and thus (1.12)^— (1.12)^ leads to an 
algebraic equation 



2(^f j A U - ^kl A ij) + ®k{jPi]l ~ $l{jPi]k ~ $j{kPl]i + $i[kPl]j + $ij{k\sPs\l] ~ ®M[i\sPs\j} = 0, 



(Al) 



and its trace gives the antisymmetric part of (2.2). By multiplying by p^ and using (2.2)-(2.5) and symmetries 
of the Weyl tensor, after lengthy calculations one arrives at an equation with no &ijki terms 



P®j[i\Pj\k] + ^jiPjiA ki + pjipji^fk + <Z>ji [3p w p\q k + P[k\iPj\i}] + 2p. 



'jl 



l[i\Pj\k] 



0. 



(A2) 



A. 2 Non-twisting mWAND with $A = o 

As already mentioned, by differentiating algebraic constraints and using the Bianchi and Ricci identities one 
can arrive at further algebraic constraints. For example, in the non-twisting case with $^ = 0, 12 by dif- 
ferentiating (3.2) we arrive at (3.4) (equivalent to setting Ay = = <3?^ in (2.3)), which can be rewritten 



as 



-[p 2 ]$ ik + Hp 2 )ik + 2<f> ij (p 2 ) jk + Qijkiip 2 )^ = (A. 



where [p 2 ] denotes the trace of p 2 , etc. The trace of (A3) vanishes and multiplying (A3) by p ik gives 



[p 2 ][* • p] - \p][* ■ p 2 } = (^ = 0=*A). 



(A3) 
(A4) 



Differentiation of (A3) leads to 



-3 



-[p 3 ]$ ife + <S>(p 3 ) ik + 2$ ij (p 3 ) jfc + * i:f -jbi(P 3 )ji 



+ ($[p 2 ] + • p 2 ]) Plk - ($p + [* • p\) (p 2 ) lk = {A*; = 0= <&£•), 



(A5) 



with trace of (A5) vanishing due to (A4). 

For = 0, the (N — l)th derivative of (2.2) reads 



®ijkl(p N )jl = °> N = l,2,... (Aij = 0= <& 



(A6) 



B Non-twisting type III Einstein spacetimes in arbitrary dimension 

For non-twisting Einstein spacetimes of type N, pij = diag(a, a, 0, . . . , 0) [8, 11], where a vanishes in the case of 
Kundt spacetimes. Let us present here a new derivation of the form of the optical matrix for a non-twisting 
type III Einstein spacetime in arbitrary dimension. This is considerably simpler than the original derivation 
of [8]. 

For non- twisting type III, eq. (59) of [8] reduces to (in the notation of [11]) 

2% jPlk = p% ij . (Bl) 

In the Kundt case (p = 0) this becomes an identity, while for non-twisting non-Kundt Einstein spacetimes 
(to which we restrict from now on) necessarily p ^ (cf. Proposition 1 of [10]) and (Bl) gives non-trivial 
information. Following [24] we choose a parallelly propagated frame with diagonal form of pij = diag(p2, P3, • • • )• 
Without loss of generality we can assume that VD^ij / f° r some values of i, j. Eq. (Bl) then implies 

P2 = \+ 0. (B2) 
12 Recall that in the non-twisting case with ^ the form of p is fully determined by point (i) of theorem 1.1, cf. [12]. 
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Using the r-dependence of the non- vanishing eigenvalues of p^ (see eq. (4.2)), this gives 



fVE^- (B3) 

p p 

where p corresponds to all non- vanishing eigenvalues of pij. Apart from the Kundt case with p^ = 0, this 
equation can be satisfied only for pij of rank 2, with a =p2 = p3- Thus the optical matrix for the non-twisting 
type III Einstein spacetime is either vanishing (Kundt) or p^ = diag(a,a,0, ... ,0), exactly as in the type N 
case. 



C On the geodetic and shearfree condition in four and higher dimensions 
C.l Goldberg-Sachs in four dimensions 

The four-dimensional Goldberg-Sachs theorem is usually expressed as a statement about a null congruence 
being geodesic (k = 0) and shearfree [a = 0). As well-known, this property can be interpreted as different but 
equivalent geometrical statements, which we now briefly review. This will later be useful for discussing the 
higher dimensional case. 

C.l.l Integrability of two-spaces 

By considering the commutator [5, D] (see, e.g., p. 77 of [3]) one immediately sees that £ is geodesic and 
shearfree if and only if the distribution 

£> = Span{m,£} (CI) 

is integrable (cf. [4]), i.e. [V,T>] C V. A similar conclusion holds of course for V. The geodesic&shearfree 
property can thus be seen as a statement about the null bivector £ Am instead of one about the vector field 
I. While the vector field £ defines a privileged real null direction lying in V n V, m is not fixed uniquely 
(null rotations with £ fixed, spins and boost can be performed, under which £ A m is fixed, up to a rescaling, 
and k = = a is preserved). The significance of this viewpoint appears clearly when considering complex 
extensions of the Goldberg-Sachs theorem [34]. 

C.l. 2 Conformal structure of the screen space 

Given a generic null vector field £, its flow preserves the conformal structure of the "screen space" L^/L 
when [35] (see also [31,36,37]) 

C e g = pg + £ (g> £ + £ ® £ o k = = a. (C2) 

In four dimensions only, this is equivalent to requiring that the complex structure of the screen space is 
preserved by the flow of £ [35-38] . 

This condition can also be derived by imposing Maxwell equations on a null field (two-form) [39] (cf. also, 
e.g., [3,4,31,36,37]). 

C.2 Conformal structure of the screen space: the standard geodesic&shearfree condition 

Ref. [35] (see also [31]) proposed to extend (C2) to higher dimensions as a definition of the geodesic&shearfree 
condition. As it turns out, 

C t g = pg + l®Z + Z®l & Ki = = S ij --^6 ij , (C3) 
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so that indeed £ is geodesic and shearfree in the standard sense. Some comments on algebraically special 
Einstein spacetimes admitting such a congruence have been already given in section 2.2 and appendix E. 

Note that for n > 4 eq. (C3) does not follow from Maxwell's equations for a null field (intended again as a 
two- form), while geodeticity (k« = 0) does, and shear is non-zero for expanding solutions [11,40]. 

C.3 Existence of an optical structure 
C.3.1 Complex notation 

For later discussions it will be useful to use a complex basis. Namely, in even dimensions one can take the 
frame (£, n, fj, A , p, A ) and in odd dimensions the frame (£, n, fj, A , p, A , x). Except for the unit spacelike vector x, 
all the vectors are null, the complex fi A are defined by 

1 1 

M 2 = ^7f( m (2) + im (3))> ^4 = ^/f( m (4) + im (5)), (C4) 

and fi A by their complex conjugates, where A = 2fi, fi = 1, . . . , (n — 2 — e)/2, with e = 0, 1 for even and odd 
dimensions, respectively (in 4d this would be the standard NP frame). The metric becomes 

g = £ (g> n + n (g> £ + fx A (g> p, A + p, A (g) fi A + ex (g> x. (C5) 

Then one can define the complex counterpart of the Ricci rotation coefficients. These are defined as an 
obvious extension of the real coefficients, e.g., 

C'AB = VA^B^bL, C AB = fl A fi b B V b l a , M.CB = Hc^B^bHAa, (C6) 

and other coefficients (and their complex-conjugates) are defined similarly. 



C.3. 2 Optical structure 

Ref. [41] studied the consequences of Maxwell's equations for a null field defined as an n/2-form in n even 
dimensions. From these, they arrived at a generalization of the geodesic&shearfree condition different from the 
k = = a condition discussed in section C.2 (see also [31,36-38]). Recently, this has been extended also to 
odd dimensions [19]. Namely, consider the totally null (n — e)/2-dimensional distribution (recall that A in fi A 
can be only even) 

V = Span{£,/x 2 ,...,/z (ra _ 2 _ e) }. (C7) 

If T> and T) 1 - are integrable, i.e., 

[V,V] C V, [V ± ,V ± ] C £> x , (C8) 

V is said to define an " optical structure" [19] (note that V 1 - = V in even dimensions). In 4d, eq. (C8) reduces to 
the standard conditions k = = a (as discussed in section C.l.l), which indeed corresponds to the conditions 
coming from the Mariot-Robinson theorem [3,4] for a Maxwell null two-form. In higher dimensions, using the 
complex counterpart of the commutators given in [9], one finds that (C8) is equivalent to 

A A A A x 

£ao = £xo = £ab = £ax = £xA = Mbo = MxO = Mbc = M.Bx = Mab = 0. (C9) 

AC A B 

(In fact one first obtains conditions such as Mbc = Mb A, but Mbc = —Mac-, etc.) The equations above 
represent the conditions obtained when n is odd, however, if n is even one can simply drop all equations 
containing x (this will be understood from now on). 

Note, in particular, that the vector £ is geodesic {Cao = £x0 = 0), but generally is not required to be 
shearfree (except when n = 4). In even dimensions it has also been observed that (C8) means that the complex 
structure of the screen space is preserved [36-38]. We have shown in [12] that in 5d a very large class of 
algebraically special Einstein spacetimes possesses an optical structure. In appendix D we extend the results 
of [12] by showing that, in fact, all algebraically special 5d Einstein spacetimes possess (at least) one optical 
structure. 
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Optical structure in six dimension Let us rewrite the conditions (C9) in real notation in the special case 
n = 6, which is useful for the discussion in the main text (see [12] for the case n = 5). One readily gets 



k 1 = (i = 2,...,5), 

P22 = P33, P23 = ~P32, /»44 = P55, P45 = ~P54, 

P24 = P35, P42 = P53, P34 = ~P25, P43 = ~P52, 
2 3 3 2 

Mm = M 50 , M 40 = -M 50 , (n = 6) (CIO) 

2323 2323 
M 4 2 - M 4 3 = M53 + M 5 2, -M52 + M53 = M43 + M42, 

223322 33 
M44 - M55 = M45 + M54, M45 + M54 = -M44 + M55- 

In particular, if £ is twistfree then p has two pairs of repeated eigenvalues. 



C.4 Optical constraint 

Based on results for Kerr-Schild spacetimes, Ref. [21] put forward yet another possible generalization of the 
shearfree condition for mWANDs in higher dimensions. This is the so called "optical constraint", already 
discussed in section 2.2, which involves only the null direction £ (as opposed to the optical structure discussed 
above). The Lorentz transformation freedom of null rotations preserving £, boosts and spins is thus retained in 
this case (indeed spins can be used to arrive at the canonical form (2.7)), see also [24] for related comments. We 
have shown in [12] that in 5d a very large class of algebraically special Einstein spacetimes admits an mWAND 
obeying the optical constraint, and Ref. [20] extended our result to prove that in fact all algebraically special 
Einstein spacetimes admit such an mWAND (see also [7] ) . Note that in 4d the optical constraint is a necessary 
condition for £ to be a repeated principal null direction but is not sufficient [12]. 



C.5 Integrability of a (complex) two-dimensional totally null distribution 
C.5.1 V23 integrable 

As a further generalization of the geodesic&shearfree condition one can consider the integrability of the complex 
two-dimensional totally null distribution 

V23 = Span{/i 2 ,£}. (Cll) 

It is easy to show that V23 is integrable if and only if 

22 2 
£20 = 0, £b2 = Mbo, £B2=Mm (B^2), C x2 = Mxo- (C12) 

In 4d this again reduces to the standard n = = a condition. 

In particular, in six dimensions (C12) can be rewritten in real notation as 

k 2 = = re 3 , P22 = P33, P23 = ~P32, 

2 3 2 3 

P42 = M40, P43 = M40, P52 = M50, P53 = M50 (n = 6). (C13) 



C.5. 2 P23 integrable with totally geodesic integral surfaces 

We can strengthen the above conditions by further requiring the integral surfaces of T>23 to be totally geodesic. 
The corresponding equations read 

£bo = = C x q, Cb2 = = C X 2, £§ 2 = (B 2), 

Mbo = = MxO, Mb2 = = M x2 , Mbo = ° (5/2), M B 2 = {B^2). (C14) 
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Now £ is necessarily geodesic (Cbo = = C x q). In 4d (C14) is equivalent to (C12) because B = 2 is the 
only possibility (and there are no x-components). In real notation (C14) can be rewritten as (where we define 
^2,3) 

«i = 0, P22 = P33, P23 = ~P32, P~ k2 = = P'k3 

M ko = = M ko , (C15) 

2 3 2 3 

M k2 = M k3 , M k3 = -M k2 . 



D Optical structures in five dimensions 

Proposition 4 of [12] gives a set of sufficient conditions for a five-dimensional Einstein spacetime of type II 
or more special to possess an optical structure [19]. In particular, it shows that, except possibly for Kundt 
spacetimes (and for a special subclass "(iii)" of genuine type II, later proven not to exist [20]), all algebraically 
special Einstein spacetimes admit an optical structure in five dimensions. It is the purpose of this appendix to 
show that this in fact holds also for Kundt spacetimes. (Obviously, if a complex optical structure is integrable 
its complex conjugate is integrable too and this will be understood in the following.) Combining this with the 
proof of [20] that the special subclass (iii) of genuine type II is empty, we arrive at 

Proposition D.l. In a five- dimensional Einstein spacetime admitting a multiple WAND £ there always exists 
an optical structure. In the case of type D spacetimes there exist in fact (at least) two optical structures. 

Proof. Let us show that Einstein spacetimes of the Kundt class always possess an optical structure in five 
dimensions. In other words, we need to show that in such spacetimes there always exists a null frame 
{£, n, m(2), wi(3), ?Ti(4)} such that the totally null distribution 

V = Span{m,( 2 ) + im^, £}, (Dl) 

and its orthogonal complement 

V L = Span{m (2 ) + im (3 ) , m (4 ) , £}, (D2) 

a 

are both integrable. This is equivalent to [12] (see also appendix C, and (1.6) for the definition of Mbc) 

Ki = 0, p 33 = p 2 2, P32 = ~P23, P2A = = P34, PA2 = = /?43, (D3) 

2 3 2 3 2 3 

M40 = = M40, M42 = M43, M43 = -M42- (D4) 

Kundt spacetimes admit a metric in the form [26,42,43] 

ds 2 = 2du [dr + H(u, r, x)du + W a (u, r, x)dx a ] + g a p(u, x)dx a dx 13 , (D5) 

where a, /3 = 2, 3, 4 in five dimensions. Here £ = d r is a geodesic, twistfree, shearfree, non-expanding mWAND, 
so that (D3) is automatically satisfied. Now, define a null frame with i a dx a = du, n a dx a = dr + H du + W a dx a 
and the spacelike vectors mu\ living in the three-dimensional transverse (Euclidean) space spanned by the x a 
(their components will be, in particular, independent of r). Then one immediately finds 

Mjo = 0, (D6) 

so that the first of (D4) is satisfied in this frame. 

Next, using the Christoffel symbols given in [26] one can check that 

m(i)a;/3 = m(i)a\\l3, (D7) 
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where the covariant derivatives on the r.h.s. is taken w.r.t. the transverse metric g a p. 
It follows that 

i i 

M jk = M jk , (D8) 

where the connection coefficients on the r.h.s. are those computed w.r.t. the transverse metric g a p (and the 
o.n. frame vectors of the transverse space are simply the "projections" fh^ a of the m^ a vectors). 
One then finds 

44 222 333 

[m (2) + im (3) ,m (4) ] = (M 24 + iM u )rh {4) + [-M42 + i(M u - M 43 )]rh {2) + [M 24 - M 42 - iM 43 ]rh {3) , (D9) 
so that Span{rri( 2 ) +im( 3 ),m( 4 )} is integrable iff 

32 32 

M 43 = M 42 , M42 = -M43. (D10) 

Proving the integrability of V and V 1 - is thus now reduced to proving the integrability of Span{m( 2 ) + 
mi(3),m(4)}. The transverse frame is arbitrary, so we just need to show that there exists at least one frame 
satisfying this integrability property. 

If we define the complex null vector field 



1 

72 v 



A* = — (™(2) +»W»(3)). ( D11 ) 



using (D10) the required integrability condition [fx, rh^] = afi + (3rh^ reads 

Ma||K 4 )M 6 = 0- (D12) 

Thus by choosing a complex null geodesies fx in the transverse three-space we automatically obtain the in- 
tegrability of the corresponding distributions V and V 1 - and (together with the result of [20]) our proof is 
complete. □ 

E Shearfree twisting spacetimes (even dimensions) 

Twisting geodesic mWANDs with zero shear are forbidden in odd dimensions [10] but they are permitted in 
even dimensions and they have necessarily det(p) / (as can be easily seen in a frame adapted to Aij, using 
the fact that Sij oc 5ij). Here we present an explicit example in six dimensions. To our knowledge, this is the 
first such example that has been identified. 

First, consider the six-dimensional Ricci flat Taub-NUT metric [44] 

ds 2 = -F(r)(dt - 2ni cos^id^i - 2n 2 cos 6» 2 d^ 2 ) 2 + 



F(r) 

+ (r 2 + n\)(&6l + sin# 2 d(/> 2 ) + (r 2 + ra 2 )^ + sin^d^), (El) 

where 

= r 4 /3 + {n\ + n 2 )r 2 - 2mr - n\n\ 
(r 2 + nf)(r 2 + n 2 ) 

We observe that this is a spacetime of type D. A geodetic mWAND is given by 

£ a dx a = dt + F(r)~ 1 dr - 2m cos Oidfa - 2n 2 cos 6» 2 d^ 2 , (E3) 
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while a second one can simply be obtained by reflecting £ as t — > —t, fa — > — fa, fa — ^ —fa [15]. Using the 
frame vectors 



m 



a dx a = ^jr 2 + n\&Qi, m^ a dx a = ^Jr 2 +n\ sin6>id0i, 



(2) 



m( 4 ) a dx a = Jr 2 + n^d8 2 , m^ a dx a = Jr 2 + n 2 . sin 9 2 dfa, 



(E4) 



one finds 



r 2 + n\ 



r 2 + n\ 



r 2 + n? 
r 

r 2 + 







n 2 



+ n 2 



\ 











n 2 



r 2 + w 2 



r 2 + 7T-2 

r 

r 2 + n| 



(E5) 
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One can easily check that p obeys the optical constraint (2.6). Moreover, we also observe that in the 
spacetime (El) the maximally totally null distribution V = Span{m( 2 )+im( 3 ) , m^+im^ , £} defines an optical 
structure (concept introduced and discussed in [19,27,31,36,37,41]), and both the totally null distributions 
^23 = Span{m,(2) + imt^,£} and P45 = Spanjm^) + imt^,£} are integrable, with totally geodesic integral 
surfaces (see appendix C for the corresponding definitions and conditions). The null vector field £ is generically 
shearing (S22 = S33 7^ 5*44 = £55), however in the special case n\ = n 2 (corresponding to the solutions of [45]) 
it becomes shearfree (while still being expanding and twisting). The class of shearfree twisting spacetimes is 
thus non-empty in higher dimensions. 

F Non-degenerate, non-twisting geodesic mWANDs violating the optical 
constraint (n > 7) 

In the main text we have seen examples of Einstein spacetimes admitting a non-twisting mWAND violating the 
optical constraint, see e.g. table 1. However, all of them have a degenerate p (i.e., m < n — 2). In this appendix 
we provide also some examples with m = n — 2, relevant to theorem 1.1 (point (hi)) and proposition 4.1 (point 
2.). 

F.l Metric 

Using the theory of conformal Einstein spaces [29] (reviewed, e.g., in [46]), by taking a double Brinkmann warp 
one can construct the following Einstein space (satisfying = 2Kg a i ) /{n — 2)): 



ds 2 = \r 2 du 2 + 2dudr + (Xur - l) 2 da 2 x + r 2 dS^, 



(Fl) 



where 



A 



A 



\n- l)(n- 2)' 



(F2) 



dc 2 is a Riemannian Einstein space of dimension n a with Ricci scalar R a = n a (n a — 1)A and dEg a Ricci-flat 
Riemannian space of dimension ns- In order for ds 2 to have a non-zero Weyl tensor, the metrics da 2 or dX 2 , 
cannot be both of constant curvature, so that either n a > 4 or rig > 4 (or both). The complete spacetime has 
thus dimension n = 2 + n a + nz>7 (unless it is of constant curvature, which is of no interest to us). 
Let us define the vector field £ with covariant and contravariant components 



£ a dx a = du, £ a d a = d r , 



(F3) 
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which is obviously null and hypersurface-orthogonal, thus automatically geodesic and twistfree, and r is an 
affine parameter along it. This can be accompanied by another null vector n (also hypersurface-orthogonal) 

n a dx a = dr + ^r 2 du, (F4) 
which satisfies the normalization condition £ ■ n = 1. 

F.2 Weyl tensor 

From the results of [29] together with those on direct product spacetimes [15] it readily follows that 

Cuabc = = Crabc- (F5) 

Let us consider a set of n — 2 spacelike o.n. vectors mu\ orthogonal to £ and n, composed of a subset of 
vectors denoted by i = A, B, . . . living in the subspace of da\ and a subset denoted by i = I, J, . . . living in the 
subspace of dEg. Then the only non-zero Weyl frame components are given by Cabcd and Cijkl, and their 
/•-dependence is (recall the first definition in (1.7)) 

&ABCD = ( Aur _ ^2 C ABCD, ^IJKL = ^C^ JKL , (F6) 

where quantities with superscript do not depend on r (these are in fact the Weyl components of the respective 
"subspaces" ) . 

In particular, one has 

= 0. (F7) 
Note also that one has generically Wij / (recall definition (3.6) and constraint (3.7)). 



F.3 Multiple WAND(s) and optics 

From the above results it follows that both £ and n are double WANDs and the spacetime is thus of type D. 
In fact a bit more than that, since 

C a bcd d = = C a bcd,n d (F8) 

(cf. [47,48] for the meaning of such conditions in terms of the Weyl type) and the type is D(abd). 13 

The optical matrix p of £ is obviously symmetric, it is moreover diagonal and non-degenerate with compo- 
nents 

Pab = t t^ab, pi j = Sij. (F9) 

\ur — 1 r 

There are thus two eigenvalue-blocks of dimension n a and ri£, or [n a , n^]. For instance, for n = 7 we can con- 
struct explicit solutions with eigenvalues {a, a, a, a, b}, for n = 8 we can have {a, a, a, a, a, b} and {a, a, a, a, b, b} 
and so on (in all cases a ^ ^ b, a ^ b). These provide examples for of the shearing spacetimes of theorem 1.1 
and proposition 4.1, although an eigenvalue structure more general than the one of metric (Fl) can in principle 
exist. In particular, the optical constraint is clearly violated (cf. eq. (2.7)). Note also that we cannot construct 
a six-dimensional example with this method. By setting the cosmological constant to zero in the above metrics 
(A = 0) one obtains Ricci-flat spacetimes - however, these are direct products and p becomes degenerate, which 
is of no interest to the present discussion. 

By taking direct products one can trivially generate many similar examples with a degenerate £. These will 
have / 0, however still with $4 = (see [15]), in agreement with point (i) of theorem 1.1. 

13 The construction of metric (Fl) is carried out more naturally in a different coordinate system (not adapted to i) in which I 
and n are related by a coordinate transformation ( "time- reflection" ) that leaves the metric invariant: therefore they share the same 
geometric properties. 
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F.3.1 On possible additional (m)WANDs 

One might wonder whether spacetimes (Fl) admit other multiple (double) WANDs (different from i and n, 
eqs. (F3) and (F4)) and in particular whether those can obey the optical constraint. Take a generic null vector 
k = ad u + f3d r + j A m^ +^ I m^ (where the m^) [ m (i)) have components only in the subspace of dcr 2 [dS 2 ,]), 
i.e., a 2 \r 2 +2a(3 +5ab7 A ^ E ' +&ul I l J = 0. Since for j A = = 7 7 k gives the directions of £ and n (respectively 
for a = and /3 = — ^a\r 2 ), we now assume that j A and j 1 are not both zero (so that a 7^ 0). Using (F5), 
(F6) and the Bel-Debever criteria [48] it follows that: (i) k is a WAND C% CD -f A j c = = C^ jkl j i j k (so 
that, using also (F5), one gets C a b c dk b k d = 0); (ii) k is an mWAND C% CD ~/ A = = C° IJKL -f I (in which 
case, using also (F5), one gets C a bcdk d = 0). 

These conditions on the Weyl tensors of da 2 and dE 2 , are generically not satisfied (see an example in 
section F.5), so in general spacetimes (Fl) do not admit any other WANDs (not even single) apart from £ 
and n. However, it is also clear that in special cases additional mWANDs may exist. For example, when 
either dcr 2 or dE 2 , are conformally flat (which is necessarily the case if n a < 4 or ns < 4) spacetimes (Fl) 
admit a continuous infinity of mWANDs (generically non-geodesic but some can be geodesic, see an example 
in section F.4). 

F.4 An explicit example with additional mWANDs 

For the sake of definiteness, one can for instance construct an explicit solution for n = 7 by taking (Fl) with 

da 2 x = V(p)dr 2 + V-\p)dp 2 + p 2 (d8 2 + sin# 2 d0 2 ), V(p) = 1 - ^ - Xp 2 , 

dSg = dz 2 , (F10) 

and a coordinate range such that V(p) > 0. 
Taking £ = d r and the orthonormal vectors 

m {2) = A- 1 V- 1 f 2 (p)d T , m (3) = A- 1 V 1 f 2 (p)d p , m (4) = A^p^de, 

m( 5 ) = A~ 1 p~ 1 sin 6~ 1 d ( p, m (6) = rl d z , A = Xur — 1, (Fll) 

one finds 

An 1 , , 

P22 = P33 = P44 = P55 = T r, P66 = ~- (f 12) 

Xur — 1 r 
The only non-zero Weyl frame components read (recall definition (3.6)) 

W 23 = W 45 = 1 -^ 4' W 2 4 = W 25 = W 3i = W 35 = -—^—-^^. (F13) 
{Xur — iy p 6 (Xur — l) z 2p 6 

However, note that since dE 2 , = dz 2 is (conformally) flat, it follows from section F.3.1 that metric (Fl) 
with (F10) admits also other mWANDs of the form k = ^[d u — \r 2 {X + j 2 )d r + jd z ], where 7 is an arbitrary 
function. 

It turns out that if 7 >r = 7 ;U = 7 z = the mWAND k is geodesic, in which case it becomes twistfree iff 7 
is a constant and it obeys the optical constraint iff 7 2 = |A| (so that it is also twistfree), i.e., 



d«-^(A + |A|)r 2 d r ±v^ 



(F14) 



In this case the corresponding optical matrix is of the form {a, a, a, a, a} for A > and {a, a, a, a, 0} for A < 
(in both cases a / 0). 
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F.5 An explicit example without additional mWANDs 

It follows from section F.3.1 that in order for metric (Fl) to admit only two double WANDs we need n a > 4 
and ws > 4, i.e., at least 10 spacetime dimensions. If we now take da\ as in (F10) but dSg corresponding to 
the Riemannian version of 4d Schwarzschild, by looking at the conditions of section F.3.1 it is easy to see that 
indeed the only mWANDS are the I and n discussed above (eqs. (F3) and (F4)). Such metric thus constitutes 
an example of an Einstein spacetime with all double WANDs violating the optical constraint. Recall also that 
such mWANDs (F3) and (F4) are geodesic and non-degenerate. 
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